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Abstract 

We illustrate the renormalized perturbation expansion method by ap- 
plying it to a single impurity Anderson model. Previously, we have shown 
that this approach gives the exact leading order results for the specific 
heat, spin and charge susceptibilities and leading order temperature de- 
pendence of the resistivity for this model in the Fermi-liquid regime, when 
carried out to second order in the renormalized interaction U . Here we 
consider the effects of higher order quasi-particle scattering and calculate 
the third order contributions to the H 3 term in the impurity magneti- 
zation for the symmetric model in a weak magnetic field H. The result 
is asymptotically exact in the weak coupling regime, and is very close to 
the exact Bethe ansatz result in the Kondo regime. We also calculate the 
quasi-particle density of states in a magnetic field, which is of interest in 
relation to recent experimental work on quantum dots. 

email address: a.hewson@ic.ac.uk 

1 Introduction 

Systems with strong local inter-electron interactions have been the focus of much 
theoretical work in recent years, as these include a variety of interesting systems 
ranging from high T c superconductors, heavy fermions and Mott insulators, to 
mesoscopic systems such as quantum dots. Conventional perturbation theory 
cannot deal with strong interactions in general, so that new techniques have to 
be developed to make predictions for the behaviour of such systems. Special- 
ized techniques, such as the Bethe ansatz or conformal field theory, have been 
successfully developed for certain classes of systems, such as one dimensional 
systems and impurity models. However, techniques are required that can be ap- 
plied more generally, particularly for systems in two and three dimensions. One 
approach which has been extended to a wider class of problems is the numeri- 
cal renormalization group approach (NRG) as developed by Wilson jjj, which 
was orginally successfully applied to models of magnetic impurities. In this ap- 
proach the higher energy excitations are progressively eliminated to deduce a 
sequence of effective models for the behaviour on lower and lower energy scales. 
The behaviour on the lowest energy scales can be calculated from the limiting 
fixed point Hamiltonian of this sequence, and its leading correction terms. In its 
original form the method only works for impurity systems but a modified form 
of the approach, the density matrix renormalization group method (DMRG) 
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^ pj, has been successfully developed for one dimensional systems. In princi- 
ple the DMNRG could be extended to two and three dimensions but there are 
technical difficulties in practice, though some calculations for two dimensional 
systems have been carried out. In another development the NRG has also been 
extended to higher dimensional lattice models by use of dynamical mean field 
theory (DMFT) §. This approach exploits the fact that certain infinite dimen- 
sional lattice models, such as the Hubbard and periodic Anderson models, can 
mapped onto effective impurity models, together with a self-consistency condi- 
tion H Q|. The calculations for the effective impurity models can be carried 
out using the NRG and iterated until the self-consistency condition is satisfied. 
The DMFT has considerably extended the potential range of application of the 
NRG approach. The NRG, however, is not the only way of realizing renormal- 
ization group ideas. The earlier way of applying the renormalization group, as 
originally developed in field theory, was via a reorganisation of the perturbation 
expansion, such that the expansion could be carried out in terms of the renor- 
malized parameters. This rearrangement of perturbation theory enabled one 
to circumvent the problem of the divergences which had plagued the conven- 
tional approaches. The elimination of the divergences, however, was essentially 
a by-product of this approach and it is possible to use the reorganisation of the 
perturbation expansion as a strategy for dealing with the low energy behaviour 
of a wide variety systems. The renormalized perturbation theory approach could 
be particularly useful in situations where there are strong renormalizations of 
the basic parameters, such as in the Fermi-liquid regime for heavy fermions, 
where the masses of the electrons may be renormalized by a factors of the order 
of 1000. 

In earlier work S we have shown how the renormalized perturbation theory 
can be applied to impurity models. In particular, we have shown that this 
approach, when applied to the Anderson impurity model and carried out to 
second order in the renormalized interaction U , gives the exact leading order 
results for the specific heat, spin and charge susceptibilities at T = 0, and 
the leading order temperature dependence of the resistivity in the Fermi-liquid 
regime. In this paper we begin with a brief review of these earlier results, and 
we then calculate some higher order corrections, in particular, we calculate the 
H 3 term for the impurity magnetization in a weak magnetic field H to third 
order in the renormalized interaction term U, and compare the results with the 
exact Bcthc ansatz result. In the final section we calculate the quasiparticle 
density of states in a magnetic field. These latter results are of some interest in 
comparing with the low temperature linear response results on quantum dots 
in the presence of a magnetic field. 



2 Renormalized Perturbation Theory 

We formulate the renormalized perturbation expansion for the impurity Ander- 
son model |9j . The Hamiltonian for this model is 

h = J2 edd ° d ° +Und ^ nd >i + J2( Vkd ° c k,<?+ v k4 : ^^ (!) 

(7 k.a k,a 



which describes an impurity d level e^, hybridized with conduction electrons of 
the host metal via a matrix element Vk, with an interaction term U between the 



electrons in the localized d state, where nd, a — d\ r d a . The width of the localized 
bound state for U — depends on the quantity A(u) = TrJ2k l^4| 2( ^( w — e k)- 
For a conduction without any prominant features this function does not have 
a strong dependence onw, so it is usual to take the case of a wide conduction 
band with a flat density of states where A(w) becomes independent of to and 
can be taken as a constant A. The partition function Z for this model can be 
expressed as a functional integral over imaginary time dependent Grassmann 
variables, corresponding to the electron creation and annihilation operators, 
using the standard methods (see for example fLOf), 

Z= I ' ]JV(d a )V(d a )V{c k ^)V{c k ^)e~ s , (2) 

J a 

where the action S is given by 

S = £am (V)Gh-, (3) 



and the Lagrangian for the Anderson model is given by 

£am(t) = 'Y^d a (T)(d T - e d )d a (T) + ^ c k , a (r)(d T - efe)c fe!CT (r) 

O k.(7 

+ f/n T (r)n i (r) +^V r fe (4(T)cfc, CT (r) + c fc)er (r)d (7 (r)), (4) 

where n CT (r) = d a (T)d a (T). One can formally integrate over the Grassmann 
variables for the conduction electrons, as they involve only quadratic terms, 
and express the result in terms of a reduced action S le a, 

Z= { \[v{d a )V(d a )e~ s ^ , (5) 

"* a 

where S rc d is given by 

S Ted = dr \ dr , V4(T)[GW(T-T0]- 1 ^(T , ) + ^ / drn^n^r), (6) 
Jo Jo a Jo 

with G^ 0) (t) = l//3X) n G ! i 0) (iw„)e- iu, ™ T , where u n = (2n + 1)tt//3. The non- 
interacting Green's function for the localized electron G&(iuj n ) given by 

G^\iLo n ) = ——- 1 (7) 

iu n + ah - e d + iAsgn(w„) 

where we have included a coupling to a magnetic field H, and h given by h = 
gfiBH/2. The Grassmann variables are required to satisfy antiperiodic boundary 
conditions d(j3) = -J(0) and d(j3) = -d(0). 

The Fourier transform of the corresponding retarded one-particle double- 
time Green's function G d ^ (cj) for the localized d electron can be deduced by 
analytical continuing to real frequencies iuj n — > Lu — i5 (5 — > +0). On introducing 
a corresponding self-energy E CT (u;, h) the interacting retarded Green's function 
can be written in the form, 

= Z h l-A V I hY ® 

U! — 6d + (711 + t/X — 2-i<j(U), h) 



In the conventional perturbation expansion this self-energy is calculated in pow- 
ers of the interaction U. It will be convenient to write the self-energy in the 
form £ CT (w + ah,h) because the non-interacting Green's functions, which are 
the propagators in the perturbation expansion, are functions of the combined 
variable to + ah. In the renormalizcd perturbation theory the perturbation ex- 
pansion is reorganized to a form appropriate for the low energy regime. The 
first step is to write the self-energy in the form, 

£ CT (w + ah, h) = £ CT (0, 0) + (to + ah)^(Q, 0) + £ r , cm (w, h), (9) 

which does nothing more than define the remainder self-energy SJj, cm (cj + ah, h), 
except that we have assumed that Luttinger's result Q that £^(0,0) is real. 
When this is substituted back into equation (||) , the Green's function takes the 
same form with a 'renormalized' energy level, width of the localized state and 
self-energy, which are denoted by a tilde, defined by 

e d = z(e d + £ CT (0,0)), A = zA, t a (u,h) = zY? a cm (to, h), (10) 

where z, the wavefunction renormalization factor, is given by z = 1/(1 — £^(0, 0)), 
and £ CT (0,0) and £^.(0,0) are to be evaluated at T = as well as to = h = 0. 
These will be the parameters of the renormalized theory instead of e d and A, 
which are specified in the 'bare' Hamiltonian of equation ([[]). Note that the g- 
factor coupling to the magnetic field H is unrenormalized. The overall z-factor 
is removed by rescaling the Grassmann fields, d a (r) — ► *Jzd a (T). 
The last parameter specifying the renormalized theory is the renormalized in- 
teraction U. This quantity is derived from the irreducible four point vertex 

function r CT ./ (to, to'), which is a special case of the more general irreducible four 

// /// 

point vertex function V a £ (to, to'; to", a/") with a" = a, a'" = a', uj" = to and 
to'" = to'. This latter quantity is derived from the two particle Green's function 
of the d-electrons in the usual way. The renormalized four point vertex function 
is defined by r CT)CT /(w,a/) = z 2 Y ^i (to, to'), and takes account of the rescaling of 
the local fermion fields. The renormalized interaction U is then defined by the 
value of T aM i (to, to') at to = to' = 0, 

C/ = f^(0,0). (11) 

As certain of the interaction effects are taken into account ab initio in the 
renormalized theory compensating terms have to be introduced to avoid over- 
counting. The Lagrangian for Anderson model can be rewritten in the form, 

C A M(a T a,d cr ,ed, A,J7) = C A u(da,da,ed, A,U) + Lcr(d a , d a , X\, A 2 , A 3 ), (12) 
in terms of the renormalized fields, where the counter-term Lagrangian is given 

by 

Ccr(d a ,d (7 , Ai, A 2 , A 3 ) = d a (T)(\ 2 d T + Xi)d a + A 3 fi T (T)n i (T) (13) 

where Ai = — z£(0,0), \ 2 = z—1 and A3 = z 2 (U — Tj j(0, 0)). By construction 
the renormalizcd self-energy £ CT (w) is such that 

£ ff (0,0)=0, £^(0,0) = 0, (14) 

so that T, a (to) = O(to 2 ) for small to, on the assumption that it is analytic at 



uo = 0. As r CTjCT (0, 0) = we also have 



rV,(o,o) = c/(i-^, CT ')- 



(15) 



The quasiparticle or renormalized Green's function takes the form 



l 



(16) 



uj — e c i + er/i + tA — S CT (a;, ft-) 



The reorganized perturbation theory is set up to calculate the renormalized 
self-energy Ecr(w, h). The propagators in this expansion correspond to the non- 
interacting quasiparticles in the Lagrangian Cam^o-, d a , 6d, A, U) with [7 = 0. 
The quasiparticle interaction U is used as an expansion parameter but all the 
terms in the counter Lagrangian Cct have to be included as well. To organize 
the expansion in powers of U the terms Ai , A2 and A3 have also to be expressed 
formally as a powers in U, 



The coefficients \\ , \\ and A 3 are then determined by the requirement that 
the three normalization conditions, (|l4|) and (|l5|), are satisfied to each order in 
the expansion. These normalization conditions are essentially those used within 
field theory in order to circumvent the problem of infinities arising from the 
lack of an ultraviolet cut-off (see for example H))- The procedure, however, 
makes no mention of infinities, it simply allows the field theoretic perturbation 
expansion to be expressed in terms of the experimentally observed masses and 
interactions. In condensed matter systems divergences do not arise in this way, 
as there is always a natural cut-off, so this is no necessity to reorganise the 
perturbation expansion. However, for the Anderson model there are very strong 
renormalizations of the effective d-level and the interactions at low energies in 
the Kondo regime, where the impurity d-electrons are virtually localized, which 
make it desirable when working in this regime to take account of these very 
strong renormalizations from the start. It allows one to develop a perturbation 
theory with an effective d-level and interactions appropriate to this energy scale. 
More generally it makes a direct link to Landau Fermi liquid theory. There are 
no cut-off dependent ultraviolet divergences to eliminate so the question of the 
renormalizability of the model does not arise. 

We can obtain significant results with this approach even at zero order, 
U = 0. If we calculate the quasiparticle occupation number nd,cr at T = and 
H = from @J) with U = 0, we find 



As the factors of z in this expression cancel it is equivalent to the exact Friedel 
sum rule [jl3| and expresses the one-to-one correspondence between the quasi- 
particle number and electron number in Landau Fermi liquid theory [ jl5| , n^ a = 
7id,cr' Hence, the d-level occupation at T = can be calculated from the zero 
order renormalized Green's function. 




(17) 




(18) 



The Freidel sum rule also holds in the presence of a magnetic field and the equiv- 
alent expression for the occupation of the d-level in terms of the renormalized 
self-energy is given by 



1 1 _ x / e d -ah + ^(0,h) \ 

n d ,a = tan = . (19) 

Z n \ A I 

We have to use the perturbation theory to calculate the field dependence of 
the renormalized self-energy. However, we can show that it is sufficient to work 
only to first order in U to obtain the exact result. There are two terms to first 
order, one from the tadpole or Hartree diagram, and one from the corresponding 
counter-term diagram in Ai. There is no wavefunction renormalization to this 
order so A2 = and also to this order ff,j,(0, 0) = U, so A3 = 0. To satisfy 
the renormalization conditions the counter-term should cancel the contribution 
from the tadpole diagram for T = h = so Ai = Un,_(0, 0). Hence the 
combined contribution is 

S«( w ,fe ! T) = i/(nf ;r (/i,r)-nf,(0 ! 0)). (20) 

The spin susceptibility of the d-electrons at T = can be calculated from 

fl , A*B(^d,T ~ n d,|)/2, by substituting the self-energy from ( |2C| ) into equation (|l9|), 
and then differentiating with respect to H . The charge susceptibility can be 
calculated in a similar way and the two results are 

where pd(0) is the quasiparticle density of states at the Fermi level and is given 

M0) = ^|j. (22) 

It is not obvious that these results to first order in U are exact. However, there 
are Ward identities |l4| ] which can be derived from charge and spin conservation, 
and in terms of the renormalized self-energy and density of states take the form, 



dh 



dfi 



= -PdA0)U. (23) 



The spin and charge susceptibilities can be derived from these exact relations 
on using ([l9|) to give 

Xd = ^j^p d (0)(l - dt/dh) = to^!p d (0)(l + ^(0)), (24) 

and 

Xd,c = 2p d (0)(l + = 2p d (0)(l - Up d (0)), (25) 

confirming that the first order results for these quantities are indeed exact. 

The impurity contribution to the low temperature specific heat coefficient 
from the non-interacting quasiparticles (U = 0) is given simply by 

Id = ^Pd(0). (26) 
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Figure 1: A plot of the renormalized parameters U and nA for the symmetric 
Anderson model in terms of the bare parameters U and nA. In the comparison 
of these parameters with 4Xk for U ^S> nA the value of 7k is given in equation 
13). 



This result corresponds to the exact result calculated by Yamada 14 . It is a 
general feature of Fermi liquid theory that the quasiparticle interactions do not 
give any corrections to the linear coefficient of specific heat as their contributions 
to the specific heat are of higher order in temperature. 

In the local moment or Kondo regime the local charge susceptibility must 
go to zero, so from equation ( [HI ) we find Upd(0) = 1. If we define the Kondo 
temperature 7k by \d — (ffMB)/47k then pd{0) = l/4Tk and all the results 
can be written in terms of Ik- They correspond to the exact results for the s-d 
or Kondo model |l6| 0. 

From the exact Bethe ansatz results 0, [TtJ for the spin and charge suscep- 
tibility for the symmetric Anderson model it is possible to deduce the renormal- 
ized parameters, A and U , in terms of the bare parameters A and U. These are 
are shown in figure 1. Initially U ~ U for small U, but when U/irA > 2, the en- 
ergy scales U and irA merge in the strong coupling regime and U = irA = 47k- 

To calculate low temperature conductivity to order T 2 one needs to calculate 
the renormalized self-energy to order cu 2 and T 2 . There is a T 2 contribution 
to the conductivity arising from the scattering of free quasiparticles by the 
impurity but there is an additional contribution due to the mutual scattering of 
the quasiparticles due to the inter-quasiparticle interactions. The lowest order 
contribution of this type arises from the second order diagram for £ shown in 
figure 2(b). The only counter term diagram that has to be taken into account 
to order lu 2 or T 2 is due to the second term in A^ 2 "* , which is required to cancel 
the linear in u> term arising from the diagram in figure 2(b), and this gives 
a!, 2 ' = 3 — 7r 2 /4. We calculate this for the case of particle-hole symmetry which 



is such that e d = U/2 and Tid(T) = 1. The contribution to the imaginary part 
of the renormalized self-energy from 2(b) is 

Imt(cj, T) = ttU 2 J p d {t)p d {e')p d (uo ~ e - e')D(u>, e, e')dede', (27) 

where 

D(u>, e, e') = (1 - /(e) - f(e'))f(e + e' - w) + /(e)/(e'), , (28) 

and /(e) is the Fermi factor 1/(1 + e^ e ) with (3 = l/T. To calculate the conduc- 
tivity to order T 2 we need to evaluate this integral to order uo 2 for T = 0. For 
T = equation ( p8| ) becomes 

D(u>, e, e') = (0(e) + 0(e'))(% - e - e')) + (1 - 0(e))(l - (9(0), (29) 

where 0(x) is the step function. To find the cu 2 coefficient we differentiate twice 
with respect to to and use the relations, 

?M = S(l) , «W_f W . (30) 

As Prf(O) = 1/nA for the case of particle-hole symmetry, the result is 

ImS(c,0) = -\U 2 PM^ 2 = (jOj . (31) 

We need the corresponding results to order T 2 for u> = 0. For = the 
temperature dependent factor in the integrand of ( p8| ) is 

D(0, e> O=2/(e)/(O(l-/(e + O). (32) 

We can change the variables of integration to x and x', where x = eT and 
x' = e'T, and the integral of equation (27) to order T 2 becomes 

rp2 / jj \ 2 r°° r°° 



Im£(0,T) = — ^- I ) / / 2F(x)F(x')(l - ,F(x + x'))dxdx', (33) 



oo •/ — oo 



where .F(x) = 1/(1 + e x ). The integration over x' can be carried out to give 
T 2 ( U V f°° x , tt 2 T 2 / U ^ 2 



lm£(0, T) = — \ — - / ^-—dx = — - • (34) 

To evaluate the conductivity we need to evaluate the transport relaxation life- 
time t(u, T) which is proportional to the inverse of the impurity density of states 
Pd(w, T) which in turn is proportional to the imaginary part of the renormalized 
Green's function so that 

/ rp\ T ,-i X t a, ^ (^-Re£(^,T)) 2 

r(w,T) cx irp d (uj,T) = A - ImEfw, T) + -= — ~- — — . 

(A — lmL(w, J )) 

(35) 



When these results are used to evaluate the contribution to the impurity con- 
ductivity <Ti mp (T) to order T 2 we find 

a imp (T) = trojl + y (T) 2 (1 + 2(R - I) 2 ) + 0(T 4 )}. (36) 

where R is the Wilson ratio given R = 1 + U/irA. This is an exact result to 
order T 2 which was first derived by Nozieres |fl8f for the Kondo regime, which 
corresponds to U = irA = 4Tk and R = 2. The more general result was derived 
by Yamada @. More recently Lesage and Saleur [jl9| have also calculated the 
coefficients of the T 4 and T 6 terms in this expansion in the Kondo regime, using 
boundary conformal field theory. 

Nothing has been omitted in the renormalized perturbation expansion, and 
it gives the asymptotically exact results in the low temperature regime, when 
taken to second order in U, so it would be interesting to extend the results by 
including higher order terms. One possibility would be to include all the terms 
to fourth order in U, and calculate the coefficient of the next correction term 
in the conductivity, the T 4 term, to compare the result with that of Lesage 
and Saleur. However, this would be require an expansion of the self-energy in 
terms of both the frequency and temperature for all the fourth order terms, 
which, though straight forward to carry out, would be a rather long and tedious 
exercise. An alternative way of examining the contributions from the next order 
terms would be to calculate the H 3 term in field dependence of the impurity 
magnetization in a weak magnetic field. The linear term in H was given exactly 
by the first order renormalized expansion. The coefficient of the H 3 term is 
known exactly for the Kondo model at T = 0, and also for the symmetric 
Anderson model, from Bethe ansatz calculations Jjl| [l7| . The renormalized 
perturbation calculation of this coefficient to order U 3 is described in the next 
section. 



3 Higher order terms 

We will perform the renormalized perturbation calculations here in a slightly 
different but equivalent way from the one used in the previous section. It will 
have the advantage of not involving the explicit use of counter-terms. We will 
also obtain an expression for the renormalized parameters in terms of the bare 
ones, at least for weak coupling. We first of all use the standard perturbation 
theory in U, and then calculate the renormalized parameters explicitly to the 
appropriate order. We can invert this relation and then write the standard per- 
turbation result in terms of the renormalized parameters, i.e. we renormalize 
the standard perturbation terms order by order using the renormalization con- 
ditions ( [l0| ) and ([ll]) . The result will correspond to the renormalized expansion 
in U, as described in the previous section, when taken to the same order. 

3.1 Third order perturbation theory 

We use the zero temperature formalism where the impurity Green's function 
can be written in the form, 

. 1 




where h = g^H/2. An expression for the impurity magnetization in terms 
of the magnetic field dependent self-energy at T = can be derived from the 
Friedel sum rule, where the impurity level occupation number nd, CT in the spin 
channel a is given by 

1 1 _ x f e d -vh + Z a (0,h) \ 
nd,a = ^ - - tan ( ^ J > ( 38 ) 

which is equivalent to equation (|l9|). We can deduce from this an expression 
for the induced impurity magnetization, and expand it to order h 3 . However, it 
will be useful to separate out the skeleton tadpole diagram shown in figure 2(a), 
which has the full Green's function, indicated by a double propagator, in the 
bubble, as this is equal to Und.-a, where rid, -a is the exact expectation value 
of the occupation number. We write the self-energy Yi a {u>, h) in the form, 

T,cr(cj,h) = Und,-<r + f}(r(w,h), (39) 



and substitute it into equation (p8|). We write the impurity magnetization 
M(h) = g^(nd^ — nd,i)/2, in a weak field as a power series, 

/ , \ 2n+l 
n ^ ' 

For the particle-hole symmetric model (e<2 = U/2) these coefficients in terms of 
the self-energy £f(0, h) is 



(1 - U/nA) V dh 



fe=0 



1 f»r3 A 2 9 3 S T (0,/i) 

A/ 3 = 771 — XT M i ' 



3(1-J7/ttA) V 2 d 3 h 



h=0 



(41) 



(42) 



We have the results for the first derivative of the self-energy with respect to h 
to order U 3 from the calculations of of Yamada H] , 



a a 



—a I —a 

4 ► 4 

a 

(a) 

Figure 3: Third order diagrams with a tadpole insertion. 

The only unknown term in the expression for the third order magnetization to 
order U 3 is the third order derivative of the self-energy at zero frequency with 
respect to the magnetic field h. 

To second order in U there is only one diagram which contributes to h), 
that shown in figure 2(b), which gives a contribution, 

zf b) (w,h) = U 2 f G\ 0) (u-Lo',h)IP^(u;',h)^;. (44) 

The particle-hole W y > ha ' propagator and the corresponding particle-particle 
propagator W y ' ha arc both evaluated in appendix A. 

where the coefficient C has been evaluated numerically, and we find a value 
C« -1.735. 

The third order diagrams fall into two types. There are three diagrams shown 
in figure 3 corresponding to dressing each of the propagators in the second order 
self-energy diagram with a simple tadpole or zero order Hartree term. For the 
symmetric model the contributions from the first two diagrams, 3(a) and 3(b), 
cancel to first order in h but contribute to higher order. The contribution from 
the diagram 3(a) in which the particle line is dressed with a tadpole is 

EfW) = v tan_1 (£) / (Gl%-^M) 2 n^V,M^. (46) 

The contribution from the corresponding diagram 3(b) in which the spin-down 
hole line is dressed is 

E< 36 >( W ,/0 = ^tan- 1 (A) / (Gf (-w W,/0) 2 n pT *V,^- (47) 

As IpTJ>i( w » = -WM(w',h) and cf } (-w + u', h) = -Gj 0) (w - u', -h), 
we can rewrite this contribution as 

= "^tan" 1 (A) J(&iH»-J-h))W M V,h)^. (48) 




We are left with the contribution from the diagram in figure 3(c) in which the 
spin | propagator of a particle-hole bubble is dressed with a Hartree bubble. 

Zf c) (u;,h) = -^tan- 1 (A) J (G f ( w - J , h)fW^{J \ h)^. (49) 

The total result from the three diagrams to order h 3 is 

^>=-K)(^) 3 +SJ© 3 . « 

where the coefficient E is calculated numerically as -5.670. 

Finally there are the two diagrams illustrated in figure 4(a) and (b), which 
can be regarded as being derived from the second order diagram, 2(b) with an 
intermediate scattering in the subdiagram corresponding to one of the dynamic 
susceptibilities. The contribution from the diagram in figure 4(a) is 

4 4( V h) = -U 3 J G?\w - «/, fc)(IP T ' fc V, h)) 2 ^ (51) 

with intermediate particle-hole scattering. For the diagram in figure 4(b) with 
intermediate particle-particle scattering the contribution is 

X ( f b) (L;,h) = -U 3 [ Gj%'-^)(n^V^)) 2 ^- (52) 
1 J 1 2m 

The total result is 

>(0, h) + £f>(0, A) = -M10 - , 2 ) {£) + g (^) 3 . (53) 

The coefficient D was estimated numerically as D=1.541. 
Collecting the results to third order in U together, 

3tt 2 277T 2 

A =16 + C*, B = 80 + C + £> + £. (55) 

4 4 

The coefficients A and B can also be deduced from the Bethe ansatz results 
for the magnetization of the symmetric Anderson model by generalizing the ap- 
proach of Horvatic and Zlatic p(J , who derived a recurrence relation for a series 
expansion in powers of U for the coefficient Mi, to obtain a similar expansion 
for M3. The details are given in Appendix B. We find A — 65/3 — 37r 2 /2 which 
gives C = 17/3 — 37r 2 /4 = —1.7355 in complete agreement with the numerical 
estimate, and B = 15(280/27 - tt 2 ) = 7.512 which agrees well the numerical 
estimate 7.515. 

3.2 Renormalization 

Having calculated all the self-energy terms to order U 3 using the standard per- 
turbation theory we want to deduce the corresponding results in the renormal- 
ized expansion to order U 3 . 





(a) 



(b) 



Figure 4: Third order diagrams with repeated particle-hole, (a), and particle- 
particle scattering, (b). 



We need to calculate the renormalized parameters, A and U, using the 
definitions given in equations ([Io|) and ( |TT|) , in terms of the bare parameters to 
order U 3 . For this we will need the wavefunction renormalization factor z and 
the irreducible four point vertex |(0, 0, 0, 0). The only contribution to z to 
order U 3 comes from the second order diagram and the result is 



U 

7rA 



(56) 



which can also be deduced from the results in Yamada's paper The con- 
tributions from the diagrams for the irreducible vertex T| |(0, 0, 0, 0) to second 
order in U cancel in the absence of a magnetic field. Diagrams which contribute 
at third order in U are shown in figure 5. There are three possible diagrams of 
the type shown in figure 5(a) each one giving a contribution U(U/irA) 2 . There 
are six other diagrams in all of the type shown in figure 5(b) and (c), and each 
of these gives a contribution (2 — ir 2 /4)U(U /ttA) 2 . The total to third order in 
U is 

1^(0,0,0,0) = U jl + (l5 - ^pj (JLy .... j . (57) 
From these two results we can deduce the renormalized parameters to order U 3 , 

...j, u = u\i-(^-^ ' 

'(58) 

These results correspond to the weak coupling region U <§; 7rA in the plot of 
the renormalized parameters shown in figure 1. 

We can invert these expression to deduce the bare parameters A and U in 
terms of the renormalized ones, A and U, 




i = A 1+ 3- 




, U = U< 1 + (tt 2 - 9) 



7rA 



(59) 




Figure 5: Third order contributions to the irreducible four- vertex IY^(0, 0, 0, 0). 



We now use these relations to express the third order result in U for the 
magnetization in terms of the renormalized parameters. It will be convenient 
to write (EG) in a modified form 



7T ^— ' V A 

n x 

The coefficient M3 to third order in U is given by 



2n+l 



M(h) = 9 -^J2M 2n+1 (i) (60) 



2 



7T 



A / \ ttA 




M 3 ={l + 4 Z )+A> Z) +B'[^ T \ }, (61) 



where 

A' = 7 + C, B' = — -4 + C + D + E. (62) 

We know from the Ward identity (|2^) that the result to order h is exact to all 
orders in U. We know that the term of order h 3 ia asymptotically exact in the 
weak coupling regime, U/ttA — > 0, so it is of interest to check it in the strong 
coupling limit against exact Bethe ansatz results for the Kondo model. We use 
the results in the previous section to express all the renormalized parameters 
in terms of the Kondo temperature 7k U /ttA — > 1 and ttA — > 4Tk as U — *■ 00. 
Wc then find 

M(h) _h_ 

5/ xb ~2T K ^192T3' ^ 

where D' = 5 + A' + B' = 8 + 5tt 2 /4 + 2C + D + £?. With the values of the 
coefficients as deduced from the third order renormalized perturbation theory 
we get D' = 12.73. If the exact Bethe ansatz result is written in the same 
form, using the same definition of Ik, the coefficient D' has the value D' = 
24V3/tt = 13.232. The error from our third order results is less than 4%. 
Hence the perturbation theory taken to third order, which is asymptotically 
exact in the weak coupling regime, is very close to the exact result at strong 
coupling. 

As U/ttA — > 1 in the strong coupling limit, the factor (U /wA) n multiplying 
the contributions from the nth order terms do not decrease with n, as they do 
in the weak coupling limit U/wA < 1. In the Appendix B we show that in the 
Kondo limit no finite set of renormalized diagrams can give the h 3 coefficient 
in the magnetization exactly. However, we have shown that the error is small 
in the limit when the perturbation series is taken to third order, and is even 
smaller for intermediate and weak coupling. This is clearly seen from the results 
in figure 6 where we plot the coefficient M3 against U/ttA over the range from 
weak (U/irA « 1) to strong coupling (U/ttA = 1) as calculated from the third 
order renormalized perturbation theory and compare it with the exact Bethe 
ansatz results, expressed in terms of the renormalized parameters U and A. 
Over the range from U = U = 0toU = 5ttA, (0 < U/irA < 0.9998), A/A 
varies by two orders of magnitude, from 1 to ~ 8 x 10~ 3 . 

It would be interesting if there would be some way of extracting the small 
h 3 correction from the higher order diagrams. We have taken account of the 
contribution from the leading order irrelevant term (U) about the low energy 
fixed point to third order, so that the remaining contributions must be related to 




U/7TA 

Figure 6: The coefficient M3 from the third order renormalized calculation 
compared with the exact Bethe ansatz result (dotted line) plotted as a function 
of U/tvA. 

the next order irrelevant terms in the effective Hamiltonian at the Wilson strong 
coupling fixed point. These terms must be a combination of local operators but 
it is not clear how to relate these explicitly to the higher order diagrams in the 
renormalized perturbation theory. 

4 Dynamic response functions in a weak mag- 
netic field 

If we calculate the w-dependence of the self-energy as well as the h-dependence 
we can deduce the form of the quasiparticle density of states in a weak magnetic 
field. If we calculate this as a general function of tu and h, rather than expanding 
in powers of u and h, it will be more convenient to revert to the renormalized 
expansion as used in section 2 with the explicit use of counter-terms. This just 
requires a rearrangement of the terms calculated in the previous section. Each 
diagram will now be interpreted as a diagram for the renormalized self-energy 
S ff (a;, h), with U — > U and A — * A, and id = —U/2 but we will have to include 
the counter-terms to order U 3 to satisfy the renormalization conditions. For the 
particle-hole symmetric model, the only non-zero counter-terms to third order 
are: A x = 0, Aj 2) = (3 - tt 2 /4), and \jp = -ttA(15 - 3vr 2 /2). The only new 
term to order U 3 is the last term which cancels off the renormalization of the the 
four vertex r^(0, 0, 0, 0) shown in figure 5, which is not needed as the vertex 
U is taken to be the fully renormalized one. An alternative way to calculate the 
counter-terms is directly from their definitions in terms of the self-energy and 
vertex functions, and to re-express these in term of the renormalized parameters. 
The counter term diagrams are shown in figure 7. The first diagram 7(a) 



•o- 

I 
I 
I 

4 

(a) (b) (c) 

Figure 7: Counter-term diagrams which contribute to the self-energy to third 
order in U. The double-dashed line represents the vertex A3, the cross represents 
the vertex A2, and the single dashed line U. 




involves the A2 vertex, and ensures that the linear term in uj is cancelled off. The 
next diagram 7(b) is an additional tadpole contribution arising from the counter- 
term interaction A3. There is also a third order counter-term diagram 7(c) 
arising from a combination of the tadpole diagram to order U with a counter- 
term vertex A 2 on the bubble. This diagram gives a contribution 



2m 



(64) 



The evaluation of the integral is straight forward and gives 

3 r 



Ef( W ,*) = A(8-2L) [Z 



tan" 1 I I 



Ah 



(65) 



If the self-energy to third order from the standard diagrams calculated in the 
previous section in terms of the renormalized parameters is written as (u>, h), 
then the renormalized self-energy to order t/ 3 is given by 



(3), 



3^ 2 



Han 




n 2 U \ ~ n 2 



h\ Ah 

tan~" I 

A h 2 



A 2 




One can check that, in the limit to = and expanded to order h 3 , this renor- 
malized self-energy, when substituted into equation (|l^), gives the same results 
for the magnetization as were obtained in the previous section. 

In figure 8 we plot the quasiparticle spectral densities in weak and strong 
coupling as a function of oj/A for h = 0.15 A. The peaks in the spectral density 
shift from w max = ±/i at weak coupling to cj max = ±4/i/3, at strong coupling. 
These are asymptotically exact results as h — * 0, as has been shown by Logan 
and Dickens . The general result for the position of the maximum in weak 
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Figure 8: The density of states fOr | and J. quasiparticles in a magnetic field 
h/A = 0.15 for weak coupling U/irA — 0.1 (dotted curves) and strong coupling 
U/ttA = 1.0 plotted as a function of oj/A. 



field from their calculation can be written in the form, 




The term in U 2 in the denominator arises from the contribution from the imag- 
inary part of the self-energy at low frequency. The peak in the spectral density 
in the strong coupling U/irA — 1, or localized limit, is the Kondo resonance, 
which has a width, A = AT^/tt. 

There have been a number of recent calculations of the Kondo resonance in 
a magnetic field for the Anderson model. The surge of interest in this topic 
has been due to the recent observations of the Kondo effect in quantum dots 
[ p2| p3[ . These mesoscopic systems can be described by an impurity Anderson 
model coupled by leads to two electron reservoirs, and tunnelling through these 
dots at very low temperatures is possible due to the presence of the Kondo 
resonance. As the parameters in the Anderson model description of a quantum 
dot depend upon gate voltages, they can be modified is a much more controlled 
way than for a real magnetic impurities, which should present a greater range 
of possibilities for comparing theory with experiment. There have been been 
numerical renormalization group calculations of the Kondo resonance in a mag- 
netic field |24|, , approximate treatments based on the Bethe ansatz equations 
[ p6[ p7{ , and also results from the local moment approach p8| . The advantage of 
the renormalized perturbation approach is that the results are asymptotically 
exact in the limits of small oj and H, and form a useful check on the other 
methods. 



5 Conclusions 



We have shown that the renormalized perturbation theory (RPT) approach pro- 
vides a way of going beyond the normal limitations of standard perturbation 
theory. The potential of this approach has been illustrated in the particular 
case of the Anderson impurity model, where we have shown that low order 
RPT calculations provide a comprehensive description of the low energy, low 
temperature range, in the Kondo as well as the weak coupling regime. Except 
in the weak coupling regime, we do not have explicit expressions for the renor- 
malized parameters in terms of the bare ones (other than than deducing them 
from the Bethe ansatz results, see figure 1), but this is also true of other ap- 
proaches to magnetic impurity problems such as the conformal field theory fl9|| . 
There may be ways of estimating the renormalized parameters using variational 
methods, numerically, or by the summing a subset of diagrams, as in the local 
moment approach which gives a good approximate interpolation from weak to 
strong coupling. The RPT approach does give a clear physical picture of the 
Fermi-liquid regime, and the low order results are asymptotically exact in this 
limit. This is also the case for other magnetic impurity models that have been 
studied, which include degeneracy and extra interactions, such as a Hund's rule 
coupling, and explicit expressions have been derived for the renormalized inter- 
actions in terms of the Kondo temperature, in the strong coupling limit p9| . 
It provides a complementary approach to the Wilson style of calculations which 



involve the explicit elimination of higher order excitations |3(J, 31. 

Two obvious questions arise: Can it be applied to lattice models? Is it appli- 
cable to systems with a non- Fermi-liquid fixed point? The method has also been 



already been extended to translationally invariant systems 32 , and related to 
Fermi-liquid theory. It is an alternative to the Wilson style of renormalization 
group approach, which has been applied to translationally invariant systems at 
one loop level by Shankar |3^] . There is potential for applications here, using 
the technique described in section 3, for applying the RPT approach to the 
strong coupling regime for such models as the Hubbard and periodic Anderson 
model. The complications that arise for translationally invariant systems and 
lattice models are mainly due to the dependence of the self-energy and renor- 
malized vertices on the wave vector k. There should be some simplification 
that one could exploit for infinite dimensional models where the self-energy is 
k-independent, and the k-dependence is suppressed at some types of vertices. 

The assumption of a finite wave-function renormalization factor z in the 
derivation of the renormalized expansion restricts our treatment to Fermi-liquids. 
Deviations from Fermi-liquid behaviour can have various causes and each case 
has to be considered on its own merits. It has proved possible to generalize 



the approach to a spinless Luttinger liquid 32 , and to the 0(3) symmetric 



Anderson model, which has a marginal Fermi-liquid fixed point |34| . 
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6 Appendix A 



There are two types of sub-diagrams which correspond to zero order spin sus- 
ceptibilities, with propagators defined by 

G^W^G^V,/^, (68) 

and 

/DO J I 

G?Hw- U ',h)GW{w',h)^. (69) 

These integrals can be evaluated analytically and the results are 

n " J "^'Ml0^) t0 ' " = ' 

A f^ii+iA^A J^±H)\ for ^o. (70) 



■kuj(uj + 2iA) { \ iA-h J \ iA + h 
and 

W^Uj, h) = -Un ( ^— ^ - — r r for oj = -2h, 

1 ' ' 2vrA \iA + h) n(iA-h) 

= £1 ^— m ( u \ iA + h ) L_ ln ( - +; A + h ) } for ^_ 

ir\uj + 2h + 2iA \ iA + h J uj + 2h \ iA-h J) 

(71) 

For particle-hole symmetry, we have W^p 1 (uj, h) = -W iM (w, h), and W l ' hi (oj, h) 
Wl> h T(w, h) = W im (uj, -h). 



7 Appendix B 

In this appendix we take the expression for the magnetization from the exact 
Bethe ansatz results of Tsvelik and Wiegmann jl7| for the symmetric Anderson 
model, and deduce a power series in U for the coefficient of the term in H 3 , 
along the same lines at that originally used by Horvatic and Zlatic ^3] for the 
order H term. The result for the magnetization to order h 3 can be written in 
the form, 



nAM(h) 



gush 

where u = U/ttA where 



( 1 -2^) e """ /!j '<"» + 2^ e3 """ /! *(">- < 72 » 



/ o poo irnrx \ 

J m (u) = ^ e-^ C -^=ld X . (73) 

" TTU J 1 - X 1 



Horvatic and Zlatic have developed a power series in u for Ji(u), 

J 1 (u) = ^G n w n , where C n = (2n - l)C n _i - ^-G„_ 2) (74) 



with Co = Ci = 1, and C2 — 3 — 7r 2 /4. The other coefficients to order u 5 are 



Q~-2 At: 2 4 _4 

C 3 = 15-^-, C 4 = 105-^ + -, C 5 = 15(63 - 7^ 2 + -). (75) 
4 4 16 16 

We can develop an expansion for Js(u) in a similar way, 

J 3 (u) = Y,C n u n where C n = (2 "~ ^ C w -i ~ X^"' 2 ' (76) 

n=0 ° 

with Co = 1, Ci = 1/3, (7 2 = 1/3 - tt 2 /4, and 

C3 ~9~Y' C '-TJ~— +16' C5 -H9 _ ~ + 16j- (7?) 
We can then write the expression for the magnetization to order h 3 in the form 

" AM C"= SC „„»-4 S ^», (78) 



5MB /i ±^ 3 A 2 



n=0 



where A„ — 3(C n +i — C„+i)/2 for n > 0. The coefficients of the terms in the 
second series to order u 4 are Aq = 1, Ai = 4, and 

65 3tt 2 . lr /280 ,\ , lr /847 91tt 2 tt 4 " 

A 2 = , A 3 = 15 7T 2 , A 4 = 15 1 

3 2 ' A V 27 / V 9 9 16, 

(79) 

In the Kondo limit the term proportional to J\{u) does not contributes to the 
h 3 coefficient, and it can be shown that the asymptotic contribution from the 
Ja(tt) term agrees with the result for the Kondo model (|63| ) with Tk defined by 

™ -7r 2 /8u+l/2u 



We can also use these results to deduce from these results the terms in 
the renormalized perturbation calculations to higher orders. We will use this 
approach to find the fourth order correction to our third order result. We can 
deduce A and U , to fourth order in U from the Bethe ansatz results for 7 and 
X- These can be inverted to calculate the bare parameters U and A in terms of 
the renormalized ones to the same order in U. The results are 

WO-K)*-^-^)*-) 

l-(l2-^)^ + 0(G 5 )). (82) 

We can then use the results above for the h 3 term magnetization to fourth 
order in U, and rewrite them in terms of the renormalized parameters. In this 
way we calculate that the correction from the fourth order terms to D', the 
H 3 coefficient in the Kondo limit, is -0.24145, which is 10% of the third order 
contribution. As the coefficient in the Kondo regime has a factor the exact 
result cannot be obtained within any finite order renormalized perturbation 
calculation, as results to finite order in U can be expressed as rational functions 
of the coefficients C n and A n , and these in turn are rational numbers and powers 
of 7r, which cannot generate to finite order the irrational number 



References 

[1] K.G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 

[2] S.R. White, Phys. Rev. Lett. 69, 2863 (1992); Phys. Rev. B 48, 10345 
(1993). 

[3] O. Sakai, Y. Shimizu and T. Kasuya, J. Phys. Soc. Jpn. 58, 3666 (1989). 

[4] Density Matrix Renormalization, ed. I. Peschel et al. (Springer, Berlin, 
1999). 

[5] A. Georges, G. Kotliar, W. Krauth and M.J. Rozenberg, Rev. Mod. Phys. 
68, 13 (1996). 

[6] Y. Shimizu and O. Sakai, J. Phys. Soc. Jpn. 65, 2632 (1996); Y. Shimizu, 
O. Sakai and S. Suzuki, J. Phys. Soc. Jpn. 67, 3181 (1998). 

[7] R. Bulla, Phys. Rev. Lett. 83, 136 (1999). 

[8] A.C. Hewson, Phys. Rev. Lett. 70, 4007 (1993). 

[9] PW. Anderson, Phys. Rev. 124, 41 (1961). 

[10] JW. Negele and H. Orland, Quantum Many-Particle Systems, Addison- 
Wesley (1988). 

[11] N.N. Bogoliubov, and D.V. Shirkov, Introduction to the Theory of Quan- 
tised Fields. (3rd edition) Wiley-Interscience (1980); L.H. Ryder, Quantum 
Field Theory. Cambridge University Press (1985). 

[12] J.M. Luttingcr, Phys. Rev., 121, 942 (1961). 

[13] J. Fricdcl, Can. J. Phys. 54, 1190 (1956); J.M. Langcr and V. Ambegaokar, 
Phys. Rev. 164, 498 (1961); D.C. Langreth, Phys. Rev. 150 516 (1966). 

[14] K. Yamada, Prog, in Theor. Phys., 53, 1970 (1975); 54, 316 (1975). 

[15] P. Nozieres, The Theory of Interacting Fermi Systems Benjamin (1964); 
A. A. Abrikosov, L.P. Gorkov and I.E. Dzyaloshinskii Methods of Quantum 
Field Theory in Statistical Physics, Dover (1975). 

[16] N. Andrei, K. Furuya and J.H. Lowcnstcin, Rev. Mod. Phys. 55, 331 (1983). 

[17] A.M. Tsvelick and P.B. Wiegmann, Adv. Phys. 32, (1983). 

[18] P. Nozieres, J. Low Temp. Phys. 17, 31 (1974): Low Temperature Physics 
Conference Proceedings LT14 5, eds. Krusius and Vuorio. p. 339. North 
Holland/Elsevier (1975). 

[19] F. Lesage and H. Saleur, Phys. Rev. Lett. 82, 4540 (1999). 

[20] B. Horvatic and V. Zlatic, J. Physique 46, 1459 (1985). 

[21] D.E. Logan and N.L. Dickens, Preprint, (2001). 

[22] D. Goldhaber-Gordon, J. Gores, M.A. Kastner, H. Shtrikman, D. Mahalu 
and U. Meirav, Phys. Rev. Lett. 81, 5225 (1998). 



[23] W.G. van der Wiel, S. De Francheschi, T. Fujisawa, J. Elzerman, S. Tarucha 
and L. Kouwenhoven, Science, 289, 2105 (2000). 

[24] T.A. Costi, Phys. Rev. Letters, 85, 1504 (2000). 

[25] W. Hofstetter, Phys. Rev. Letters, 85, 1508 (2000). 

[26] J.E. Moore, X.-G. Wen, Phys. Rev. Letters, 85, 1722 (2000). 

[27] R.M. Konik, H. Saleur and A. Ludwig, \cond-mat/010304& (2001). 

[28] D.E. Logan, M.P. Eastwood and M.A. Tusch, J. Phys. Condens. -Matter, 
10, 2673 (1998): D.E. Logan and N.L. Dickens, Europhysics Letters, 54, 
227 (2001), 

[29] A.C. Hewson, The Kondo Problem to Heavy Fermions (Cambridge Univ. 
Press, Cambridge, 1993 and 1997): J. Phys. Cond. Mat, 5, 6277 (1993). 

[30] O. Sakai, Y. Shimizu and T. Kasuya, J. Phys. Soc. Jap. 58, 3666 (1989). 

[31] T. A. Costi, A.C. Hewson and V. Zlatic, J. Phys. Cond. Mat. 6, 2519 
(1994). 

[32] A.C. Hewson, Advances in Physics 43, 543 (1994). 
[33] R. Shankar, Rev. Mod. Phys., 66, 129 (1994). 

[34] R. Bulla, G-M Zhang and A.C. Hewson, Phys. Rev. B, 56, 11721 (1997): 
S.C. Bradley, R. Bulla, A.C. Hewson and G-M Zhang, Euro. Phys. Journ. 
B, 11, 535 (1999). 



